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IN THE CLAIMS - Following is the list of claims and their status: 

1 . (Original) A method for determining from measured reflection data 
on a plurality of trace positions one or more subsurface parameters, said method comprising the 
steps of: 

(a) preprocessing the measured reflection data into a plurality of partial or full 

stacks; 

(b) specifying one or more initial subsurface parameters defining an initial 
subsurface model; 

(c) specifying a wavelet or wavelet field for each of the partial or full stacks 
of the measured reflection data; 

(d) calculating synthetic reflection data based on the specified wavelets and 
the initial subsurface parameters; 

(e) optimizing an objective function, comprising the weighted difference 
between measured reflection data and synthetic reflection data for a plurality of trace positions 
simultaneously; and 

(f) outputting the optimized one or more subsurface parameters. 



2. (Original) The method according to claim 1, wherein the step of 
optimizing the objective function comprises minimizing the objective function: 



# stacks # traces 



F residual ~ ^ ^ residual t i ^ Lp ires idual{Sij— Mij) 
»=1 j=\ 
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wherein sy represents trace j of measured reflection data for stack i, my represents trace j of the 
modeled synthetic reflection data for stack i, w reS iduai,i represents a weighting factor for stack i, 
#traces represents the total number of traces, #stacks represents the total number of stacks and 
Lp,residuai is an adjustable norm of the residuals (sy - my), which provide the deviation of the 
synthetic reflection data from the measured reflection data. 

3. (Original) The method according to claim 2, wherein the objective 
function comprises one or more stabilization terms and/or one or more correction terms. 

4. (Original) The method according to claim 3, wherein a stabilization 
term is a measure for the deviation of the reflectivity away from 0. 

5. (Original) The method according to claim 4, wherein said measure 

comprises 

Us tacks # trace s 

F reflectivity = ^ reflectivity, i L ,p \ reflectivity (ji.j) 

wherein ry is the reflectivity trace for stack i and trace j, w re fi e ctivity,i is a weighting factor for stack 
i, #traces is the total number of traces, #stacks is the total number of stacks and Lp, re nectivity is an 
adjustable norm of the reflectivities. 
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6. (Original) The method according to claim 3, wherein a stabilization 
term is a measure for the parameter contrast. 



7. (Original) The method according to claim 6, wherein said measure 

comprises 



U pa rame ters tt trace s 
F contrast ~~ / . ^ contrast, n Lp contrast (Cj.n) 



wherein Cj, n is the contrast trace for the n th subsurface parameter, w cont rast,n is a weighting factor 
for the n th parameter, #traces is the total number of traces, #parameters is the number of 
parameters, and Lp jCon trast is an adjustable norm of the contrasts. 



8. (Original) The method according to claim 3, wherein a stabilization 
term is a measure for the deviation of the subsurface parameters from the initial subsurface 
parameters. 



9. (Original) The method according to claim 8, wherein said measure 

comprises 



# parameters 



F initial ~~ ^ W ' initial, n Lp initial (jPj,n~ P initial J,n) 
n=\ 7 = 1 
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wherein (pj >n - Pinitiaij,n) represents the trace with the difference between the calculated subsurface 
parameter n at trace position j and the corresponding initial subsurface parameter, Wi n i t i a i,n is a 
weighting factor for the n th parameter, #traces is the total number of traces, #parameters is the 
number of parameters and L P ,i n itiai is an adjustable norm of said difference. 

10. (Original) The method according to claim 7, wherein a stabilization 
term is a measure for the deviation of the calculated subsurface parameters from a priori specified 
functional relationships between subsurface parameters. 

1 1 . (Original) The method according to claim 1 0, wherein said measure 

comprises 

# f uncti ons U trace s 

F functions = ^functions, v ^ Lpj unc tionsf v (Pj, . . . JPj, # parameters) 

v=l y=l 

wherein f v represents the deviations of the subsurface parameters at trace j away from the v th 
functional relationship between different subsurface parameters, Wfunctions^ is a weighting function 
for the v th functional relationship, #traces is the number of traces, #functions is the number of 
functional relations and Lp^nctions is an adjustable norm of said deviations. 

12. (Original) The method according to claim 3, wherein a stabilization 
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term is a measure for the lateral variability of the parameters. 

13. (Original) The method according to claim 12, wherein said measure 

comprises 

# parameters # n eighb ors U trace s 
Fiaterat= ^ S S ™ lateral* Ojj) Lpj atera t (djj n ) 

n=\ 1=1 j=\ 

wherein dj jn ,i is the difference of the samples of parameter p n at traces j and 1, corrected with any 
difference in the initial model, wi at erai,n (Tjj) is a trace for parameter n describing the weighting for 
each parameter sample where this weighting is a function of x^i which is a trace which at each 
parameter sample provides a measure of the local correlation between the traces j and 1, #traces 
is the number of traces, #neighbors is the number of neighboring traces used in the calculation, 
#parameters is the number of parameters and Lp,i a terai is the adjustable norm of said differences 
<W 

14. (Original) The method according to claim 13, wherein the parameter 
difference djj, n is defined as 

(djXn){tk) = PlMk + Atj.l.k) ' Pj.nitk) - ip initial, l t nifk + A//,/,*) ~ P initial jM id) 

wherein Atjj jk is the time shift at parameter sample k which time aligns the parameters of trace 
1 to trace j at sample k, where surrounding trace samples are interpolated if at time tk + Atj^i a 
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sample is not defined, where 13,1 is now the local correlation incorporating the time shift and 
Lp,iaterai is an adjustable norm on the parameter differences. 



15. (Original) The method according to claim 1, wherein a correction 
term is a measure for the differential time shifts between traces of measured reflection data 
stacks. 

16. (Original) The method according to claim 15, wherein said measure 

comprises 

U stack s n trace s 

Ftime = ^ ^ Lpjimei? ij " T gij) 

1=2 j=\ 

wherein To,ij is the trace with the initial time values of the time stretch and squeeze control points 
for stack i and trace j and T id is the time of shifted control points, w t i me ,i is a weighting factor for 
stack i, #stacks is the number of stacks, #traces is the number of traces and L P ,time is an adjustable 
normalization factor of the difference between x 0 ,ij and Tjj. 



1 7. (Original) The method according to claim 1 , wherein a stabilization 
term is a measure for the parameter differences between reflection data acquisition surveys taken 
at different points in time. 
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18. 



(Original) 



The method according to claim 17, wherein said measure 



comprises 




# pa rame ters # trace s 




wherein (Pj, n> k-pj,n,k-i) is the difference in parameters of trace j and parameter n between survey k and 
its time preceding survey k-1, Wparametens,!! is a weighting factor for each parameter n, w survey ,k is a 
weighting factor for each survey k, L P)t imeiapse is an adjustable norm of said difference, #surveys 
represents the number of surveys and #parameters represents the number of parameters. 



and/or constraint functions are applied to one or more of the subsurface parameters; and/or are 
applied to control changes of subsurface parameters between surveys taken on different times; 
and/or constrain outside a specified subsurface zone the parameter changes between surveys to 
small values relative to changes expected within the specified subsurface zone; and/or constrain 
the minimization by setting the subsurface parameters outside a specified subsurface zone from 
survey to survey at the same value. 

20. (Original) The method according to claim 1 , wherein for optimizing 
the objective function outside a specified subsurface zone only one or one set of different 
subsurface parameters are specified. 
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The method according to claim 1, wherein constraints 
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21. (Original) The method according to claim 1, comprising the 
generation of quality control information. 

22. (Original) The method according to claim 21, wherein the quality 
information includes at least one of: (i) synthetic data based on the optimized subsurface 
parameters; (ii) residual data obtained by subtracting the synthetic data from the measured 
reflection data; (iii) deviation data obtained by determining the deviations away from the initial 
subsurface parameters; (iv) deviation data obtained by determining the deviations away from the 
corresponding functional relations; and (v) deviation data obtained by determining the deviations 
away from well log data. 

23. (Original) The method according to claim 1, wherein the reflection 
data is at least one of seismic data and time lapse data. 

24. (Original) The method according to claim 1 , wherein the subsurface 
parameters comprise elastic parameters. 

25. (Original) The method according to claim 24, wherein the elastic 
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parameters comprise pressure wave velocities and/or shear wave velocities and/or densities in 
the subsurface and/or the subsurface parameters comprise any mathematical relation between 
pressure wave velocities and/or shear wave velocities, and/or densities. 

26. (Original) The method according to claim 1, wherein the subsurface 
parameters comprise compositional parameters representing the rock and fluid composition of the 
subsurface. 

27. (Original) The method according to claim 25, wherein the seismic 
data comprises at least two seismic partial or full stacks containing different angle dependant 
information on seismic reflections in the subsurface. 

28. (Original) The method according to claim 23, wherein the time lapse 
data at each survey time comprises at least one seismic partial or full stack. 

29. (Original) The method according to claim 2, wherein the adjustable 
norm L P of a variable x for variable x comprises 

# s amp les 

S NT 
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wherein P is a user adjustable factor. 

30. (Original) The method according to claim 29, wherein the norm L P 
is normalized by or with at least one of: (i) exponentiation with 1/P; (ii) the number of 
samples; and (iii) the square root of the variance of x. 

31. (Original) The method according to claim 23, wherein the seismic 
reflection data is determined from at least one of the following source-receiver combinations: 

P-wave source and P-wave receiver, P-wave source and S-wave receiver, S-wave 
source and P-wave receiver, S-wave source and S-wave receiver. 

32. (Original) The method according to claim 1, wherein the reflection 
data is echo-acoustic data and the subsurface is human or mammal tissue or any other material. 

33. (Original) A method for determining from measured reflection data 
on one or more trace positions a plurality of subsurface parameters, said method comprising the 
steps of: 

(a) preprocessing the measured reflection data into a plurality of partial or full 

stacks; 
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(b) specifying a plurality of initial subsurface parameters defining an initial 
subsurface model; 

(c) specifying a wavelet or wavelet field for each of the partial or full stacks 
of the measured reflection data; 

(d) calculating synthetic reflection data based on the specified wavelets and 
the initial subsurface parameters; 

(e) simultaneously optimizing an objective function, comprising the weighted 
difference between measured reflection data and synthetic reflection data, for said plurality of 
subsurface parameters and for each trace position separately; and 

(g) outputting said plurality of optimized subsurface parameters. 

34. (Original) The method according to claim 33, wherein the step of 
optimizing the objective function comprises minimizing the objective function: 

# stack s # trace s 
F residual = ^ ^residual.i ^ Lp, residual (Sij—rHij) 
i=l j=\ 

wherein Sij represents trace j of measured reflection data for stack i, my represents trace j of the 
modeled synthetic reflection data for stack i, w res iduai,i represents a weighting factor for stack i, 
#traces represents the total number of traces, #stacks represents the total number of stacks and 
Lp,residuai is an adjustable norm of the residuals (sij - my), which provide the deviation of the 
synthetic reflection data from the measured reflection data. 
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35. (Original) The method according to claim 33, wherein the objective 
function comprises one or more stabilization terms and/or one or more correction terms. 



36. (Original) The method according to claim 35, wherein a stabilization 
term is a measure for the deviation of the reflectivity away from 0. 



37. (Original) The method according to claim 36, wherein said measure 

comprises 



tt stack s tt trace s 
F reflectivity = ]^ ^ reflectivity r i L <p [reflectivity (?i,j) 
i=l j=l 



wherein r;j is the reflectivity trace for stack i and trace j, w re fi e ctivity,i is a weighting factor for stack 
i, #traces is the total number of traces, #stacks is the total number of stacks and L P>re fi e ctivity is an 
adjustable norm of the reflectivities. 



38. (Original) The method according to claim 35, wherein a stabilization 
term is a measure for the parameter contrast. 



39. (Original) The method according to claim 38, wherein said measure 

comprises 
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# pa rame ters ft trace s 

Fcontrast = ^ Wcontrast,n ^ Lp, contrast! (Cj n ) 
n=l j=\ 

wherein Cj, n is the contrast trace for the n th subsurface parameter, w contr ast,n is a weighting factor 
for the n th parameter, #traces is the total number of traces, #parameters is the number of 
parameters, and Lp >CO ntrast is an adjustable norm of the contrasts. 

40. (Original) The method according to claim 35, wherein a stabilization 
term is a measure for the deviation of the subsurface parameters from the initial subsurface 
parameters. 

41 . (Original) The method according to claim 40, wherein said measure 

comprises 

# pa rame ters # trace s 

F initial = Winitial.n ^ L p initial (jP ] j ,n~ P initial J, n) 

wherein (pj, n - Pinitiaij,n) represents the trace with the difference between the calculated subsurface 
parameter n at trace position j and the corresponding initial subsurface parameter, Wi n i t i a i, n is a 
weighting factor for the n th parameter, #traces is the total number of traces, #parameters is the 
number of parameters and Lp, in itiai is an adjustable norm of said difference. 
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42. (Original) The method according to claim 39, wherein a stabilization 
term is a measure for the deviation of the calculated subsurface parameters from a priori specified 
functional relationships between subsurface parameters. 



43. (Original) The method according to claim 42, wherein said measure 

comprises 

# f uncti ons U trace s 
Ffunctions = W ? functions ,v Lpj unct i ons fj^j^\ y .,.^j^p arame ters) 

v=i y=i 

wherein f v represents the deviations of the subsurface parameters at trace j away from the v th 
functional relationship between different subsurface parameters, Wfunctions.v is a weighting function 
for the v th functional relationship, #traces is the number of traces, #functions is the number of 
functional relations and Lp,f unc tions is an adjustable norm of said deviations. 



44. (Original) The method according to claim 35, wherein a stabilization 
term is a measure for the lateral variability of the parameters. 



45. (Original) The method according to claim 44, wherein said measure 

comprises 
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# pa rame ters # neighbors U trace s 

Flateral = J] S S W/tf/era/,* Qjj) Lpj atera l (djXn) 

n=\ l=\ y=l 

wherein dj jn ,i is the difference of the samples of parameter p n at traces j and 1, corrected with any 
difference in the initial model, wi ate rai,n (xj,i) is a trace for parameter n describing the weighting for 
each parameter sample where this weighting is a function of x h \ which is a trace which at each 
parameter sample provides a measure of the local correlation between the traces j and 1, #traces 
is the number of traces, #neighbors is the number of neighboring traces used in the calculation, 
#parameters is the number of parameters and Lp,i a terai is the adjustable norm of said differences 

46. (Original) The method according to claim 45, wherein the parameter 
difference dj,i, n is defined as 

(dj,l,n)(tk) = Pl.nih + &tjj t k) ~ Pj.nifk) - (p initial Xnifk + A //,/,*) ~ P initial J,n(tk)) 

wherein Atj^k is the time shift at parameter sample k which time aligns the parameters of trace 
1 to trace j at sample k, where surrounding trace samples are interpolated if at time tk + Atj,k,i a 
sample is not defined, where rjj is now the local correlation incorporating the time shift and 
Lp,iaterai is an adjustable norm on the parameter differences. 

47. (Original) The method according to claim 33, wherein a correction 
term is a measure for the differential time shifts between traces of measured reflection data 
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stacks. 



48. (Original) The method according to claim 47, wherein said measure 

comprises 

U s tack s ft trace s 

Ftime = ^ W t i me j ^ Lpjime ( T ij ~ T Q ij) 
i=2 j=\ 

wherein x 0 ,ij is the trace with the initial time values of the time stretch and squeeze control points 
for stack i and trace j and Tjj is the time of shifted control points, Wti me ,i is a weighting factor for 
stack i, #stacks is the number of stacks, #traces is the number of traces and Lp >t i me is an adjustable 
normalization factor of the difference between x 0 ,ij and xy. 

49. (Original) The method according to claim 33, wherein a stabilization 
term is a measure for the parameter differences between reflection data acquisition surveys taken 
at different points in time. 



50. (Original) The method according to claim 49, wherein said measure 

comprises 

Us urve ys # pa rame ters # trace s 

Ftimelapse = ^ f survey \k ^parameters.n Lptimelapse (Pj.n.k- Pj.n,k-l) 

k=2 n=\ j=\ 
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wherein (pj, n ,k-Pj,n,k-i) is the difference in parameters of trace j and parameter n between survey 
k and its time preceding survey k-1, w param eters^ is a weighting factor for each parameter n, w sur vey,k 
is a weighting factor for each survey k, L P , t imeiapse is an adjustable norm of said difference, 
#surveys represents the number of surveys and #parameters represents the number of parameters. 

51. (Original) The method according to claim 33, wherein constraints 
and/or constraint functions are applied to one or more of the subsurface parameters; and/or are 
applied to control changes of subsurface parameters between surveys taken on different times; 
and/or constrain outside a specified subsurface zone the parameter changes between surveys to 
small values relative to changes expected within the specified subsurface zone; and/or constrain 
the minimization by setting the subsurface parameters outside a specified subsurface zone from 
survey to survey at the same value. 

52. (Original) The method according to claim 33, wherein for optimizing 
the objective function outside a specified subsurface zone only one or one set of different 
subsurface parameters are specified. 

53. (Original) The method according to claim 33, comprising the 
generation of quality control information. 
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54. (Original) The method according to claim 53, wherein the quality 
information includes at least one of: (i) synthetic data based on the optimized subsurface 
parameters; (ii) residual data obtained by subtracting the synthetic data from the measured 
reflection data; (iii) deviation data obtained by determining the deviations away from the initial 
subsurface parameters; (iv) deviation data obtained by determining the deviations away from the 
corresponding functional relations; and (v) deviation data obtained by determining the deviations 
away from well log data. 

55. (Original) The method according to claim 33, wherein the reflection 
data is at least one of seismic data and time lapse data. 

56. (Original) The method according to claim 33, wherein the subsurface 
parameters comprise elastic parameters. 

57. (Original) The method according to claim 56, wherein the elastic 
parameters comprise pressure wave velocities and/or shear wave velocities and/or densities in 
the subsurface and/or the subsurface parameters comprise any mathematical relation between 
pressure wave velocities and/or shear wave velocities and/or densities. 

58. (Original) The method according to claim 33, wherein the subsurface 
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parameters comprise compositional parameters representing the rock and fluid composition of the 



subsurface. 



59. (Original) The method according to claim 57, wherein the seismic 
data comprises at least two seismic partial or full stacks containing different angle dependant 
information on seismic reflections in the subsurface. 



60. (Original) The method according to claim 55, wherein the time lapse 
data at each survey time comprises at least one seismic partial or full stack. 



6 1 . (Original) The method according to claim 34, wherein the adjustable 
norm L P of a variable x for variable x comprises 



LM = 



Usamples 



wherein P is a user adjustable factor. 



62. (Original) The method according to claim 61, wherein the norm Lp 
is normalized by or with at least one of: (i) exponentiation with 1/P; (ii) the number of samples; 
and (iii) the square root of the variance of x. 
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63. (Original) The method according to claim 35, wherein the seismic 
reflection data is determined from at least one of the following source-receiver combinations: 

P-wave source and P-wave receiver, P-wave source and S-wave receiver, S-wave 
source and P-wave receiver, S-wave source and S-wave receiver. 

64. (Original) The method according to claim 33, wherein the reflection 
data is echo-acoustic data and the subsurface is human or mammal tissue or any other material. 

65. (Currently Amended) A device for determining from measured 
reflection data on a plurality of trace positions one or more subsurface parameters, the device 
comprising: 

(a) input means for inputting at least the measured reflection data and one or 
more initial subsurface parameters defining an initial subsurface model; 

(b) processing means for: 

(i) preprocessing the measured reflection data into a plurality of partial 

or full stacks; 

(ii) specifying a wavelet or wavelet or wavelet field for each of the 
partial or full stacks of the measured reflection data; 

(iii) calculating synthetic reflection data based on the specified wavelets 
or wavelet fields and the initial subsurface parameters; and 

Page 21 of 30 

{W01 48850.1} 



Application No. 10/672,751 

Paper Dated: November 4, 2004 

In Reply to USPTO Correspondence of June 4, 2004 

Attorney Docket No. 702-03 1290 

(iv) (a) optimizing an objective function, comprising the weighted 
difference between measured reflection data and synthetic reflection data for a plurality 
of trace positions simultaneously; or 

(b) optimizing an objective function, comprising the weighted 
difference between measured reflection data and synthetic reflection data for a plurality 
of trace positions simultaneously and one or more stagilization stabilization terms and/or 
one or more correction terms; and 

(c) output means for outputting optimized one or more subsurface parameters. 



66. (Currently Amended) A device for determining from measured 
reflection data on one or more trace positions a plurality of subsurface parameters, said device 
comprising: 

(a) input means for inputting at least the measured reflection data and one or 
more initial subsurface parameters defining an initial subsurface model; 

(b) processing means for: 

(i) preprocessing the measured reflection data into a plurality of partial 

or full stacks; 

(ii) specifying a wavelet or wavelet field for each of the partial or full 
stacks of the measured reflection data; 

(iii) calculating synthetic reflection data based on the specified wavelets 
or wavelet fields and the initial subsurface parameters; and 

(iv) (a) simultaneously optimizing an objective function, comprising the 
weighted difference between measured reflection data and synthetic reflection data, for said 
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plurality of subsurface parameters and for each trace position separately; and or 

(b) simultaneously optimizing an objective function, comprising the 
weighted difference between measured reflection data and synthetic reflection data for a plurality 
of subsurface parameters and one or more stabilization terms and/or one or more correction 
terms; and 

(c) output means for outputting optimized subsurface parameters. 

67. (NEW) A method for determining from measured reflection data 

on a plurality of trace positions one or more subsurface parameters, said method comprising the 
steps of: 

(a) preprocessing the measured reflection data into a plurality of partial or full 

stacks; 

(b) specifying one or more initial subsurface parameters defining an initial 
subsurface model; 

(c) specifying a wavelet or wavelet field for each of the partial or full stacks 
of the measured reflection data; 

(d) calculating synthetic reflection data based on the specified wavelets and 
the initial subsurface parameters; 

(e) optimizing an objective function, the objective function comprising the 
weighted difference between measured reflection data and synthetic reflection data, and one or 
more stabilization terms and one or more corrections terms, for a plurality of trace positions 
simultaneously; and 

(f) outputting the optimized one or more subsurface parameters. 
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68. (NEW) A device for determining from measured reflection data on a 
plurality of trace positions one or more subsurface parameters, said device comprising: 

(a) input means for inputting at least the measured reflection data and one or 
more initial subsurface parameters defining an initial subsurface model; 

(b) processing means for: 

(i) preprocessing the measured reflection data into a plurality of partial 

or full stacks; 

(ii) specifying a wavelet or wavelet field for each of the partial or full 
stacks of the measured reflection data; 

(iii) calculating synthetic reflection data based on the specified wavelets 
or wavelet fields and the initial subsurface parameters; and 

(iv) optimizing an objective function, comprising the weighted 
difference between measured reflection data and synthetic reflection data, and one or more 
stabilization terms and/or one or more correction terms for a plurality of trace positions 
simultaneously; and 

(c) output means for outputting optimized one or more subsurface parameters. 



{W0148850.1} 



Page 24 of 30 



